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ABSTRACT 


ON THE NATURAL TRANSVERSE MODES OF A SUBMARINE 


by 
EDWARD JOHN OTTH 
and 


PETER WRIGHT WOOD 


SUBMITTED TO THE DEPARTMENT OF NAVAL ARCHITECTURE 
AND MARINE ENGINEERING ON MAY 23, 1955, IN PARTIAL 
FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF 
NAVAL ENGINEER. 


A graphical method of determining normalized curvature 
patterns, deflection patterns, and frequencies for the natural 
modes of interest of a ship treated as a non-uniform, free-free, 
transversely vibrating beam is presented. Corrections for shear 
deflection and rotatory inertia are integrated into the method, 
which is based on procedures developed by Stodola and 
Schadlofsky. Integrations are performed mechanically by an 
integraph. The method is demonstrated by determining the first 
three horizontal modes of a three-eighths scale model of a 
submerged submarine, 


The method is found to be convergent, yielding acceptable 
results, and taking account of the secondary effects considered. 
Three cycles of iteration for each of the first three non-trivial 
natural modes of the model yielded frequencies of 417, 838, and 
1391 cycles per minute, respectively. Patterns obtained are tabu- 
lated in detail. 


The process developed is recommended for consideration by 
small research facilities with limited budgets. 


Thesis Supervisor: STEPHEN HARRY CRANDALL 
Ph.D. 
Title: Associate Professor of Mechanical Engineering 
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NOMENCLATURE 
total cross-section area of hull structure considered effective 
in bending and shear, (sq. inches) 


axis translation of assumed or calculated deflection for orthe- 


gonality with rigid bedy translation 


rotation coefficient ior axis of assumed or calculated deflection 


for orthogenality with sigid Lody rotation 


adjusiment coefficlenta ior correcting assumed or calculated 


deflections fer orthogonality with previous modes 
constant of integration apglied to shear (equal to sem) 
constent of integration applied to bending (equal to zero) 
constant of integration applied to slope 

constant of iategration applied to deflection 

Young's «modulus, 30 x 19 sai 

shear modulus, 11.5 x 19° psi 


moment of inertia of huli structure about vertical axis through 


centroid of each sectian, 


ratio of average shear stress to maximum shear stress over 


a cross-section for horizontal transverse shear. 
length overall of ship, ft. 
bending moment, ft-lb, 
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NOMENCLATURE (Continued) 


LX ~ 
(f - x) ) ax 
an operational constant, = [f i “es a § . ve 


shear force, lb 
period of vibration, milliseconds 
t 
an operational constent = Sv m 9 dx /L 
© 
T fie. 
an operstinnal constant = rh ¥, mg ax/ L 


deflection patterns for a uniform prismatic bar obtained from 
consideration of simple bending theory only, according to Lord 
Rayleigh. 


frequency of vibration, cps 


force of gravity, 32.2 ft/sec” 


virtual mass per unit length, including structure, ballast, fme- 
flooding water, etc., plus water entrained in horizontal vibration 


(ib see’ ft") 
m Qs virtual weight per unit length (lb/ft) 


coordinate parallel to longitudinal axis of ship, measured from 
bew (it) 


coordinate of longitudinal centroid of virtual mass 


coordinate parallel to beam of ship, measured from vertical 


plane through centroid of ship; measure of horizontal deflection 


the deflection pattern obtained from a cycle, vefore orthogona - 
lization 
value of unnormalized deflection curve taken at the normalizing 


point 
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NOMENCLATURE Continued 


=z coordinate parallel to vertical axis of ship measured from 


plane through centroid of ship. 
= virtual displacement (pounds); see 


=z component of slope due to shear deflection 

= component of slope due to bending 

= component of slope due to bending before axis determinatim, ¥+G 2-0 
= deflection before determination of axis, Hte,F Y 

= slope before determination of axis; 7 + sf 

= dumriny variable, replacing x 


= circular frequency, (radians/second) 


” 


Note: Subscript to deflection indicates mode number. 


Hyphenated subscript indicates mode number and iterative cycle num- 
ber, 
Prime indicates differentiation with respect to x. 
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I INTRODUCTION 


A. Histerical 

For the last seventy years, the subject of vibrations of ship hulls 
has attracted ever increasing attention from the naval architect, Larly 
investigators, notably Otto Schlick,emphasized the empirical treatment; 
but so complex is thearucture of a large ship, and so varied its design 
features, that it became clear that this approach would not iully serve the 
need. Ultimately then, we muat look to the theoretical approach for the 
Geeshapenew of adequate formulations. 

Basically, the ship is considered a slender free-f{ree elastic beam, 
vibrating transversely to ite leng axis, for which the theory stems [rom 
Lerd Rayleigh. [6]* From energy considerations he developed the basic 
equation based on simple Lending. He then showed that for various 
assumed mode shapes, the one for which the energy balance gave the lowest 
frequency was the correct shape. He also published a moditied equation, 
accounting for the effect of rotatory inertia, but, considering thia ef.ect 
negligible, did not emphasize it, 

Timoshenko formulated the shear flexisility term for Kayleigh's 
basic equation, and this set down the well known Timoshenko mechanisin, 
which considers the contribution of bending, shear, and rotatory inertia, 
While this for:nulation is the most complete available, it is complex and 
difficult to solve exactly, even for a uniform beam, andin recent years a 
number of papers have veen devoted to this solution, generaliy involving 
Laplace transforms or series solutions. 11, 33, 34, 63, et.al. 

Interest inthe non-uniform Leam inspired work along difterent lines, 
since a beam with a general distribution of mase, inertis, and section area 
is not amenable to exact solution. Cme approach is to attempt to jiorssulate 


one or more of these distributions inte expressions suitable for analytic 


* Numbers in brackets refer to references listed in Appendix F. 
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treatment by operational calculus, e.g. (9, 10, 35, 37] . However, it 

is evident that a reasonably accurate method of solution which is unre- 
stricted in its application, i.e., can treat completely general distributions 
of structural characteristics, has advantages over any isolated analytical | 
solution, however rigorous. Thus we shall consider the iterative methods, 
both graphical aad tabular. 

Stedola's work on turbines [8] suggested treating non-uniform bearne 
by integrating Nayleigh's difierential equation numerically or graphically. 
His methods also inciuded the use of the funicular polygon, and Kiyleigh's 
rule of minimum frequency, In using the successive integration methed, 
he observed that while the two-noded mode behaved well, a three-noded 
starting assumption diverged during iteration until it returned to the two- 
noded mode; he failed to recognize the mutual orthogonality of sormal 
medes. Viauello introduced this method into general use in the field of 

engineering, and J. =. Taylor and Schadloisky applied it to ships (illustrated 
in [66] , as is {ewis's tawular methed), All worked primarily with the 
first, or two-noded mode, and employed the simple-bending forim of the 
equation of Kayleigh, though iaylor recognized the importance of | imeo- 
shenko's shear term in ship vibrations. Koch [2] proved the convergence 
of the method, and more important, dernonustrated the necessity of 
rigorously maintaining orthogonality with all previous modes throughout 
the iterative method. 

Murray demonstrated the importance of transverse shear in shib 
vibrations [70], and introduced this effect into Schadlofsky's method, 
(rather than applying it as a subsequent correction to frequency as did 
MeGoldrick [68] )» .This refined method, first employed by Murray and 
the authors in work on a submerged submariaze [71], was noted also for 


a method of determining the slope axis by integration by parts of the 
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expression for orthogonality with rigid bedy rotation (pitch and yaw), and 
for calculation of higher modes than the first, facilitated by application of 
Koch's intermode orthogonality conditions. 

Other graphical methods have been employed by Southwell, Gumbel, 
Morrow, and others, but were suitavle generally for the first mode only, 

Recently the principles employed in these iterative methods have veen 
applied to lump pararneter formulations for solution by sequence controlled 
digitel calculators or by electrical analogy. [56, 67, 69]. In particular, 
Prohi and Myklestad developed finite difference equations for use in the 
jaylor Model Basin digital computer, taking account of simple bending, 
shear deflection, and rotatory iaertis, Calculations made b, this method 
for vibration in the vertical plane for the attack transport U.S.5. NIAGARA 
with and without the secondary effects (shear flexibility end rotatory inertia) 
were compared with experimental data. The increasing importance of 
secondary effects with increasing frequency was demonstrated, but for this 


vessel rotatory inertia was relatively unimvortant in the vertical plane, 


B., Objec tive 


This investigation is concerned specifically with transverse vibration 
in the horizontal vlane of a submerged submarine, Kecent interest in the 
bodily response of a submerine to dynamic loading such as that accompanying 
@ non-contact underwater explosion has. led to the development of metheds 
of modal analysis. These methods, and other demands of modern shi? 
design, have underscored the inadequacy of older established methods of 
determining natural niodes of response of a ship. 

In particular, then, the object of this investigation was to establish 
and demonstrate a method of determining, as accurately as possible, the 


curvature pattern, deflection pattera, and natural {requencies of the normal 
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nodes of interest in a transversely vibrating submarine of general mass, 
inertia, and section area distribution. Further, the imethod was to ve 
capatle of execution wilhout resart to such expensive mechanisms as 


digital computers, 


C. Method: 


The method arrived at is basically similar to that of Schadlofsky 

[68], involving iterative cycles operating on an assumed deflection 
pattern, each cycle comprising the four successive integmtions suggested 
by Rayleigh'’s equation. The Timoshenko shear deflection correction is 
made to the siove ovtemed from bending, and the rotatory inertia correc- 
tion of Kayleigh is made to the derivative of the moment. Koch's ortho- 
gonality checks are rnade to filter out components of all lower modes, and 
Murray's slope axis determination is employed. Integration is performed 
by integraph (filate I), and calculations by commercial desk calculator. 
The method is demonstrated for the first three normal modes, using 
structural data from an existing three-eighths scale model of & submarine, 
including an allowance for entrained water, 

Shear lag, and the possiole interaction setween shear deflection 
and rotatory inertia are not considered. 

It is noted t at the problem is formulated as an eigenvalue prab- 
lem, in which solutions exist only for a discrete set of squared frequencies, 
each of which corresponds to a particular mode pattern, except that the 
zero eigenvalue indicates both transistional and rotational vibration of 


the rigid body. 
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WT] PROCEDURE 
The basic theory of vibration of a non-uniform, transversely 


vibrating beam as formulated by Lord Kvyleigh, considering only the 


simple vending effect, can be expressed briefly by: 


_ 3? ( ard} ) = - myo" (1) 
om ox’ 


where the asaumption of simple harmonic motion of circular frequency 
has been incorporated, Four successive integrations with respect to 
this equation comprise the framework upon which the iterative rnethed of 
this investigation has been developed, ‘the deflection and frequency 
included in'the inertia loading term on the right side of (1) are assumed 
values, and the deflectioa indicated on the left, obtained by the quadruple 
integration, is the first estimate of the true deflection pattern. ‘ihe first 
estimate of the square of the true circular frequency is obiained Ly 
comparison of the assumed and derived deilection antinodal amplitudes, i.e., 


( o&* ) ( w*y*) start (2) 


4 
—_— (y*) final 


Note also that the second integra! of inertia loadiag, which is the curve of 

bending moment, must be divided by the curve o. iocul flexural rigadity 

before the last two integrations are performed, since Klis a variable in x. 
We first modiiy the process by applying to the shear curve yielded by 


the firat integration a correction for rotatory inertia, as follows: 





dM = Q - Imy (3) 
dx A 


where 7 is the bending component of slope. Since this quantity occurs 
later in ecny iterative cycle than M, we must use the 7 value from the 


previous cycle as a best approximation, ‘Thus this correction is not made 
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until the second cycle of each mode, 

The second modification of the process is the correction of the slope 
curve obtained by integrating the curvature distribution. This correction 
introduces the effect of shear flexibility, and the total slope resulting is 


expressed by: 


yee +f (4) 


where 


P= Q (5) 


Here the shear vahz © has been obteined vefaore the slope calculation, 
and the correction can be made in each cycle, 

Application of boundary conditions and orthogonality conditions, 
including orthogonality (with respect to mass) with the zero eigenvalue 
motions and with all modes pelow the one being obtained, at the 
appropriate points in and after each cycle, determine the four constants 
of integration, and refine the pattern obtained from the iterative cycle. 

In this favestigation, since mass and flexural and shear rigidities 
have general, hence irregular distributions, analytical integration was 
not possible. The process was therefore carried out in each instance by 
mens of a mechanical integraph with a 23" carriage beam (manufactured 
in Switzerland by G. Coradij. This inatrument draws an integral curve 
when its stylus is made to trace the argument curve (Plate I). Ali caleu- 
lations have been carried owt and plotted at intervals ot oue-~fortieth of 
ship length, as this was found to be the greatest luterval which would 
detect importani variations in the functions used. 

The sequence of operations is asiollows: 

(1) The structural characteristics needed are computed fron. 


available plans and tabulated (Appendix E), 
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(2) A tabulation of functions of structural characteristics and of 


(3) 


(4) 


(5) 


(6) 


(7) 


the centroid of virtual mass which will be used throughout 
the work is prepared, (Table 1V, Appendix C). To this 
table are added a tabulation ‘of a function of each modal 
deflection pattern ovtained, for usein orthogonaliziag sub- 
sequent modes, 
A tabulation of constants computed for repeated use is prevared 
(Table V, Avpendix C). 
A tabulation of the Rayleigh mode patterns (obtained from 
application of his simple-bending theory for a uniform slender 
beam) is prepared, for all modes to be obtained, These are 
normalized patterns. 
The Rayleigh pattern for the mode cveing calculated is orthogona- 
lized with translation, rotation, and all previous modes (there 
would be none for the two-neded mode), In practice, the Rayleigh 
pattern need not be used, any curve af the correct general snape 
is suitable ior a starting guess, and the better the guess, the more 
rapid the conversion to an acceptable solution. Perhaps the best 
choice would be the patiera fron. a similar ship, it such is 
available, 

The resuliing orthogonalized pattera will now te used as a 
starting point of the iterative process, 
Inertia loading based on this pattern is computed and plotigd as a 
function of X/L along the length of the ship. 
This is integrated to obtain shear, to which is added the rotetury 
inertia term, giving the derivative of bending moment, (In the 
first cycle of each mode shear ia the diifereatial of bending 


morment), M' is plotted. 
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(3) M' is integrated, and the integral curve is divided, at each 


station, by flexural rigidity to obtain curvature, which is 


plotted, 


(9) Curvature is integrated to obtain the bending component of slope 


to within a constant. 


(10) Shear is divided by shear rigidity, and the resulting shear slope 


(11) 


(12) 


is added te the vending co.cponent, and the total is used ia an 
auxiliacy process to determine a constant which locates the 

axig (the constant of the integration in (9)), which is added also, 
to give total slepe . This is plotted, 

Slope is integrated to give the deflection pattern to within an 
additive constant, which is used in an auxiliary process to obtain 
the axis shift. The resultant deflection pattern is tabulated, 
This pattern is now submitted to the same orthogonality checks 
used on the starting guess, The corrected curve is normalized 
to give a pattern whose value is unity at the center - or center- 
most -after-antinode, The unnormalized value at this antiaode 
is retained for use in the frequency calculation, The normalized 
pattern becomes the sasis of the inertia louding ior the next 
cycle, or is aceepted as the final pattern, Vhis decision is based 
on the arneunt of change of detiectioan pattern and frequency be- 


tween this cycie and the previous one, 


(13) Frequency is obtained by equation (2), where ( y*)stary = 1, and 


(1) 


(15) 


(y*) final is the unnormalized antinode deflection retained in (12). 
In the final cycle of each mode, the curvature pattera is 
normalized ty dividing throughout by (y*) fimel and retained as 
the final curvature pattern. 

The shear siope distribution is subtracted from the corrected 


total slope distribution to give the corrected bending slope. hia 
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is normalized by division by (y*) final and retained for 
use in the rotatory inertia correction of the subsequent 
cycle, 

(16) The pattern obteined in (12) is now ueed to repeat the 


entire cycle, if deemed necessary. 


Appendix A, Part 1 contains the derivations of the important relations 
used, Part Z is a detailed outline of the procedure, containing formulations 
of all operations as used in the application of the method illustrated in 
Appendices B aad C. 

ihe method was executed using etendvasel data for a three-eighths 
scale model of 4 submarine, and three lateral modes were computed, each 
having been carried through three cycles of iterations, this being sufficient 


to show convergence oi the rnethod and to reduce residual error, as 


indicated roughly by the last increment of change, to an acceptable level. 
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Hl KESULTS 


iy osservation of the convergence of the resuits of iteration, 
pattern?, normalized to give 2 value of unity at the ceater or 
centermost after antinode, are presented in Figure XXXII. The 
 orresponding curvature patteras, froai which strain patterns can 
be obtained, were normatiszed in such @ riauner that when the factor 
wetween norinalized deflection and true deflection imaxinova has been 
obtained, this aanse factor multipiled by the normalized curvature pattern 
will give the maxima of the true curvatuyve sattern, all di:nensioas 
being consistent. These patterns are presented in Figure XAXZIV. 

Table I gives the final frequencies and periods of the three modes 
caiculated for the model}: 
FABLE | 


NATURA. TRANSVERSE FREQUENCIES AND PERIODS 
FoR i09.5 FOOT SUDMEKGLD SUGMARINE MODEL 


Mode Natural F req. (cpr) Ratio Period (milliseconds) 
i 4i7 i 143.9 
2 &34 | 2,61 Ti. 
= 1391 3. 34 43.3 


* See supplementary discussion, Appendix D 

Tabie Il iists the offsets of the normalized deflection and curvature 
patterns pictured in Figures AXXUHI and XXAIV. 

Ag to the method developed, it was found to be convergent, though 
it appears that the delayed introduction of rotatory inertia inte the pro 
éess moust actto delay convergence. This is compensated for by the 
{frequent orthoyonality checks, especially their application to the 
Starting pattera. 


The results and method wre discussed in the subsequent section. 
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IV. DISCUSSION OF RESULTS 


S. Mode Patterns: 


The shape of the mode patterns (Figure XXXIIL) gemerally met the 
autho:s' expe tations. In ships and submarines, voth the masse and 
fiexu.al rigidity are greater at the midsection than at the ends. This 
has the effect of aitering the mode shape so that the center deflections 
acre smaller relative to the end deflections than in a uniforsm beam, 
and the nodal points are closer to the :nidship section, For example, 
for the first mode, the centexs deflection maxiroum of the uniform: bean: 
amounts to 60% of the end deflection, while in the case of the suv marisze 
studied, the anti-node deflection was 29. 4% of the bow deflection and 
22.8% of the stern deflection. The nodes for the uniform beam occur 
at 22% frorn the ends, while for the subtnarine investigated, they 
occurred at 26.2% from the bow and 30.7% from the stern respectively. 

Figure I clearly demonstrates the proname increase in mass, 
sectional area, moment of inertia, and shear factor K in the vicinity of 
the hard tank and conning tower amidships in the case of the submarine. 
One can readily arrive at two consequences of this discontinuity of 
structural strength, The first of these has already been mentioned - 
that of the very high end deflections relative to those in the middle vody. 
This behavior will obviously occur for all the modes, with the resulting 
trend towards patterns with many small loops in the :niddle with large 
end branches. It is easy to see how any loading, such as an explosion, 
occurring in the middle body of a subrmarine can incur very large defiec- 
tions at the ends with attendant small deflection near the load. In the 
case of the third mode, the deflection at the stern wiil be six times that 


atthe center! The reader wil) readily appreciate the importance * 


this characteristic of the mode patterns with respect to the bodiiy 





aS - 


etuvesd to HUlesooRe Vi 


reacene abe» a 


wd! joe Chietowng (HLA KA ovop! UL aavetiog be of} to eqeds ott 
bos eke 245 docu eedlimrdve baw eqiae of eewlled eg e® ‘te oadine 
eetT ubos sit 04 xual wolrsschin adi Us yetheny Sta qiiblgty isyeewsl 
sennseflad vaswes 4) Sai) oe eqads ahow adi quiliverhe bo soothe oat and 
need crvaling & ot eeds Geustyelie® bev aa) Ht eritale” atilwme et 
sivmane vot qultvee qitatin ad) G2 apcely ee Belltoy ladon a: Raw 
seed wireline od © eisuninem owlloditeh veteet 243 eben serM 04) vel 
seiemepe way lo 00s edt al bil celselie® hus pai to HPS ot sievome 
bois welibelisd wae Sat le HO 25 caw anttvelind eben vas oH lire basta 
sycoe mased ori eilae afi vol esdon ad7 eetteetla’ opens wat ln 4 58 
qadi Dategivecrat sarieedes 3d] (6b siiew chan wal) qray @E5 be 
isp lieeqaey awnse oti wou) AT DE hae wey BA} most AS .46 4 Bewsaene 

pest at Seeetani foerusnmne edd shintionnersd yisabls } wvvg!t 
le Wiaisiy oft ot & soled) Vande Boe .aitvewt fo lemmas sete lanciine® 
entudendive 0di Yo SE 00 Ot egiddblalh whwet gilnons Bae does Bred #02 
to Yrlwatteonsih lal ld eh amnepanane Set be owlvte yilbase eae SO 
. bewolsaec, weed Yhasile end sebdt bh soni act dtgasvis lan uiowsie 


joDah egral que susal B62 adthendan © le yhod wihiie edt #! grirsonse 
su ol ete) nat thew nolsveltad Sams tombaetts uti chew wits 54 eal 
cod? womeed cla ov (he nye pdt 50 weliveliod eas wees bey ots Ww fae 


= woaieanowe! 248 staluerqge ylibées tie wahery a47 | wehees of) 38 
ive db ch tanger Mitw anaeney simon sil tee silelvarceseds she 





«R65 


response of submarines. Severe whipging nay be expected to occur 
whenever the contricutions of the odd and/or even nuinvered mode end 
deflections become additive for any general vibratory motion of the 
submarine | 

The second consequence of the irregular distrivution of structural 
characteristics becomes manifest in the effect of rotatory ineatia and 
shear flexibility. Fora given structure mass per unit length, a sub- 
inarine wili be maay times stiffer than a solid prismatic bar, since 
in the case of the former, the crossection is essentially a circular ring 
ef material relatively distant frorn its ceatroid, This fact, coupled with 
that of the nature of the iongitudina! districutien of fiexural rigidity, 
explains why the magnitude of the rotatory inertia correction to the 
frequency of the first mode is sormewhat greater for the subsmarine than 
for the unifor:n beam or even a surface ship. its increased moment of 
inertia for a given mass, ft is compared to the uniform beam results in 
larger values of im . ven for equal slopes, the rotatory inertia 
correction would , be larger for the suvruarine. However, the in- 
creased fiexibility of the end sections over their uniform beam counter- 
parts causes the slove to increase over ali the length of the submarine, 
the greatest increase occurring between the ends and the nodes, The 
steepness of all the deflection patterns near the ends results ina rather 
selective effect of rotatory inertia upon the patterns, Im generai, the 
effect of rotatory inertia upon the mode patterns studied was most pro- 


nounced at the ends, and resulted in an increase in the relative defiection.- 


A comparison of experimental and theoretical mode patterns for a similar 
submarine where rotatory inertia was neglected in the latter celculation, 
shows larger experimental end deflections than the theory predicted. 


The authors feel that such a discrepancy may be due to the application 
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of an incomplete theory, particulary in the light of the results of this 
thesis, 

Shear fiexibvility apparently causes little change in the shape and 
magnitudes of the nor:nalized mode patterns. The shape of the KAG 
curve (Figure I) indicates that the location of maximium shear near 
the quarter points of the subrnarine will produce the largest value of 
shear slope, Actually the percent contribution to y’ of shear slose 
was much rnore uniform: over the full length of the submarine than was 
expected. This accounts for the small effect of shear deflection on 
the normalized patterns. As was expected, the value of @ increases 
with increasing modes nun.bers, and thus the effect of shear deflection 
upon frequency becomes more important with increasing mode number, 
Frequency tehavior is discussed more fully in a subsequent seciion, 

It should be mentioned, however, that with the third mode, the percen- 
tage of shear component of slope gave indications of becoming larger 
in the region of the middie body, where the greatest effect on ww’ will 


result. 


B, Curvature Patterns: 


The curvature patterns, Figure AAXIV, clearly indicate possible 
regions of structural weakness with respect to resistance to damage from 
Whipping. In particular, it appears that excessive strain can be expected 
atavout X¥/i. = 0.9 in the event of third mode whipping of significaat 
magnitude*® The magnitude of the strain wire, of course, also depends 
upon the radius of the submarine hull at this location. 

* Note that the rudder, propellers end shafting add mass with little 


attendant increase in stiffness in this part of the sul marine. 
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Except for the apparent deficiency just noted, the structural design 
of the submarine appears to te quite well talanced. The general shape 
of the curvature patterns obtained py the authors agrees quite favorably 


with experimentaily determined curvature patterns of similar submarines. 


CG. Frequencies: 

Frequency vcehavior observed in this investigation generally confirmed 
expectations based on theoretical consideration of the secondary elastic 
effects, Exact breakdown of these effects is not available in that no com- 
plete iterative caiculation was conducted without these effects, Although 
the first cycle of each mode involved no rotatory inertia effect, it cannot 
be considered as yielding a useful figure for comparisons, tCecause of the 
predominating effect of the starting estimate curve. However, if we accept 
these first cycle frequencies for a qualitative examination (Tabdie IL, 
Appendix I), we see the expected drop in frequency when rotatory inertia 
is introduced, of magnitude greater than what can be attributed to con- 
vergence alone. It further appears to be of increasing importance with 
higher modes. 

Data showing the independent effect of shear flexitility is not available, 
but preliminary investigations su:nmarized in refereace [71] showed a 
Sirnilar depression in frequency when this effect was introduced, although 
the effect was separated for the firat mode only. This is generally .on- 
sidered to be of increasing importance for higher modes. 

The cumulative result of these effects is noted in the fact that exper - 
irnentally determined natural frequencies correspond more nearly to the 
series |,2,3, than to the 32, 52, 7°, series obtained theoretically from 
a uniform bearn without secondary effects. [69] This same tendency is 


ovserved in the frequencies obtained in this investigation, whkh were seen 
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in Table 1 to correspond to the series |, 2.61, 3.34, or, for contrast, 
expressed in the second forin; f. 4.25%, 5. 2°. 

However, comparison with results of earlier computations for the 
same rnodel, but without the rotatory inertia effect, cast some doubt on 
the accuracy of the third mode frequency of 1391 cpm, indicating it is 
perhaps too high. This suspicion is reinforced by the 3.34 value in the 
above series. The error (probably amounting to less than 4%) is con- 
sidered to ve due to the fact that the intermode orthogonality check as 
formulated is inexact for other than the trivial modes, and the error 
becomes of some importance for the third and higher modes, The more 


correct formulation is developed in Appendix D, and the third mode 


frequency corrected by its partial apstlication. 
D. Sources of &rror: 


Several sources of error arise from fields the detailed consideration 
of which is beyond the scope of this thesis, Sut are noted for completeness. 

The first of these stem frornm the information included as Structurai 
Data, found in Appendix ©. A suvmarine, though more regular in design 
than a surface ship, is a complex structure, and the participation of its 
coraponents in carrying the loading are more difficult to analyze than 
simple beam theory suggests. Therefore agreement is hard to find as 
to which memoers houid be inciuded in calculating moment of inertia, 
cross sé€ctional area resisting shear, and the shear distrivution coeffi- 
cent K. 

Far more controversial is the matter of entrained water; hezce it is 
especialiy difficult to verify a given theory because vibration theory is 
not sufficiently refined to permit separating out other effects which inter 
act with the manifestations of entrained water. 

Hoth of the above questions have been sidestepped by utilizing data 


already calculated from plans, based on the assumptions outlined in 
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Appendix E. It is not to be inferred that these matters are insignificant; 
they are vital to any prediction of vibration, and are the subject of 
censideravle study. 

It has seen assumed that deflections are small and that the pro- 
portional limits have not teen exceeded. Standard vaiues of Young's 
Modulus and Poisson's Ratio have been assumed, Any departure from: 
these assuit:ptions ina ship tested experimentaliy can be expected to 
cause discrepancies, 

The second group of errors can be traced to the procedure and 
to ita execution. 

First, failure to carry the iteration through sufficient cycies leads 
teanerror, The residuai crror after a given number of cycles wiil be 
large if (a) the starting assumption of deflection was grossly in error, 
or (vb) if orthogonality checks were not properly and regularly applied. 
An @¢rror was introduced aiso by the application of the inexact orthogono- 
nality formulation. 

Next, the rotatory inertia term, veing based necessarily on data 
from a previous cycle, is siightly in error, 

Further, the failure to make this correction in the first cycle re- 
sults in a discrepancy which will take several cycles to eliminate entirely. 

Shear lag, and probabiy other lesser effects, have not ceen 
considered. 

Integraph operation invoives some small error. This is periodi- 
cally filtered out Sy coundary conditions and other checks, but some 
residue necessarily remains. In particular, a sinall amount of lost 
motion was noted in the integrach used. 

Any error in developing a given mode pattern will result in a smal) 
error in each subsequent pattern, through the intermode orthogonality 
checks, which necessarily employ the ordinates of previously determined 


modes, 
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Finsily, the usual remarks concerning human error in computing 
vaiues, and in plotting, integrating, and reading curves. It is noted, 
however, that the method by nature tends to reject such errorsor to 


make them immedtately ovious as the technique becomes familiar. 


E. Applicabilit : 

The cnethod of this paper is not envizgayed as replacing the rather 
advanced techniques of digital computers and electricai analogs, It 
may well complement them to some degree, but is interided primarily 
to serve the needs of the emaller research facility which, operating on a 
lirmited budget, is faced with the necessity of performing such determinations 
with simpler tools, Integraphs of the type used are found in a number of 
educational inétitutions, design offices, and shipyards, tut the authors 
have reason to velieve thet more than a few of therm lie neglected beneath 
at least a decade's accumulation of dust. The only other tool needed is an 
automatic or serni-automatic commercial desk calculator, It is considered 
that the work can be set up on appropriate forms in such a manner that 
a mathematics aid or junior research assistant can carry it out without 


special training and ander only a modicum of supervision. 
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Be 
Ve. CGCONCLUSBIONS 


Convergence of the iterative method, with secondary efiects, 
as developed in this thesis appears satisfactory, and any tendency 
towards instability is held in check by periodic application of orthogona- 
lity conditions, The method inherently detects and filters out errors, 
so that a mistake in an early cycle does not invalidate subsequent cycles 
of the same mode, out may delay convergence slightly. 

While the procedure appears complex, it can be programmed 
on standard forms for application by research assistants without special 
knowledge of vivrations or advanced mathematics. kquipment required 
is standard to many installations and inexpensive relative to such special 
equipment as digital cormputers. 

Results of this method appear to be acceptatie; the third mode 
frequency is considered accurate to within 2, 6% after three iterative 
cycles, while earlier modes, having attained better convergence in the 
same number of modes, are considerably more accurate, 

The secondary efiects are considered to have seen adequately 
accounted for, with little added complication to the process, in particu- 
lar, rotatory inertia appears to have had 4 sornewhat greater effect than 
has been evidenced in studies of uniform beams; end deflections are seen 
to Se increased relative to antinodal excursions when the eiiect is con- 
sidered. Secondary etfects thus are included not only in irequency 
determination, sut also in deflection and curvature patterns. 

& clear indication of highly losded sections is presented by the 
normalized curvature distributions which result directly {rom the pro- 
cess. Structural inadequacies are thus pointed eur. 

It thus appears that the method developed is well suited to the 
purposes of small facilities conducting vibration research on limited bud- 
gets, where accurate theoretical determination of information on natural 


modes is needed. 
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33 
VI RECOMMENDATIONS 


It is recommended that the method for modal calculations which 
is described in this thesis be considered sy small research activities 
for similar work. As a general rule, such activities operate under a 
rather limited budget with correspondingly small allocations for corm- 
puting equipment of the analog or digital types. On the other hand, most 
of these activities will certainly own more than one desk calculator and 
will have the funds available forthe purchase of an integraph in the event 
that one is not owned. At least for the modes of int¢rest, the method used 
uy the authors appears rnore than satisfactory from the standpoint of 
moderate equipment costs, simplicity, and availacility. Digital compu- 
ters are in ever increasing de:mand, and the delay involved in obtaining 
results from a central calculating contractor may be unacceptable for 
the smmali research organization. The method is clearly not competitive 
with that of the digital cornputer when the latter is readily availatbie and 
a large number of mode patterns is desired. 

&xperimental corroboration of the authors' results would be 
highly desirable before the method is accepted for general application. 
For any given submarine design this can ve accomplished by mechanical 
excitation of a rnode) of the proposed vessel, srovided it is of sufficiently 
large scale to assure structural sirmilerity. 

lt is further recommended that the effect of rotatory inertia as 
well as shear deflection be included in the natural node calculations 


regardless of the :..ethod used, 


In any future application of this method, it ia recommended that 
the more exact orthogonality conditions which arise from the inclusion 
of rotatory inertia effects be applied. While it is more convenient to 


negiect them, such practice may result in the lose of refinements, espec - 


jally in the higher modes, which are of the sarne order of magnitude as 


L40Li one cJ , 


doisw smolialeclas Labow set bedtem 42) gals beleewqmnces a! 3 

eeWivives Jocaecet (ise beembindwe Be eleads eid et Bedivoaed tt 
a wake Sinteqe Be/heas doue ,elut Lorem # Oe eeu valhcnie ool 
nqves so) sackianolie \imewe yiyalbamanes 1s > ster beghod Sasimdt sectior 
jaro ved Vedte 8) > wove! (atigib to gelame of) we ixeccqdape 9targ 
pew voeulenlae USD Rae Taes e2G00 wee Yialerees Lim webieese pend bo 
ureve adi of dqargptok oe ter sesdtonce attend? pitalione sacl edtawed lie 
bees badinsn aft JeEEMod 16 eeGoes ad) U6) mami te bewee sam 41 ene ped? 
: hs Selonbents 04) aren? yiniadis see cailt prom wtampys eAmaaie oe? od 
eqs Laatg ee AUlitealievs ace .qialiqper eters aereiapt oprsemont 
 pptesasan ad devlenat yeink afr kc Demued gniverrach ave ah sseeue! 
ot silebyeodaty B¢ Yau Tel2ensseo petteinctes Jeniibe & aren’ silanes 
eyiittepees ad (inane 6) Buceum ed7 oe lino ina, 2s COLaeke, lemme echt 
bee féalings Vibert sf teal aap nade amteqees Lediglh oti dest dite 
12) Ot ee  (heeetende at aogetiog Sia me Tekmen age 6 

od Bluew etioees ‘evedies a4) 10 seDevEset=es Lege 
emtanslinggs Lameeny ot beeg eons «t bediaw fy eecteks olven leek yidyid 
Umytamdobe: ceo bedtatlquresas we ane O20? <giees au amuve aovtg yoo ot 
qitesiatlies to a) tf tebivere leeeee beeeqove #4) iy Leta « bo endlanoes 
; wt pi elionle levees ie oeeee af sle2e sgTel 

=o slsvont -redeser te sankle-ail adh Sohecerewsns segmabenis 
gmotsatevies aber beusian att of babuisal ws en tpresiem secdece jee 
dae wrommrve: ee hee tt @ eo olira ages 

tant teleemerc cant OF 1) Dadian obit te apie pepe re tet gre al 

corwsbanh eth mews RAIA Cree Remtions Cilerreedlis Inet wer a 
oy besuseraee fume of at ote _Dellgue Oo pinnae aumead rotwicn be 
Ge? grmaransins a aa0t oe? oh luce, sam SHE here Aone wnods brelyse 
ae sharing do cebea eune est te ne maddy conker srtald O61 oF glial 


’ = 


34 


the rotatory inertia corrections themselves, 
The effect of shear lag in ship-type structures upon their vibration 


behavior should te investigated. 
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APPENDIX & 


PART 1 


Theoretical Basis: 


The subrnarine is intially considered to behave as a iong, thin 
vearm with vanishing spear force and bending moment, Unader this assump- 
tion, then, the ciassicai theory of the lateral vibration of a free-f{ree non- 
uniform bar applies. [6] [9] . We ere concerned here with the normal 
modes of vibration, When a var vibrates in one of its nermal modes, the 
deflection at any iocation varies harmonicaily with time. If damping is 


neglected, the lateral displacement of the bar can be represented by ; 
y = & ( A cos w€ + 8B sin wt ) (6) 


where x is a function of the cocurdinate x, and determines the ghape or 
pettern of the norms! or natural node of vit:ration which we are considering. 
X is known in the literature as a "normal" or natural” function . 

It is obvious that any particular normal mode af vivration is charac- 
terised soiely by x; and the frequency of that mode, oa; . Gy super- 
imposing all possible normal modes of vibration, the general expression 
for the free lateral vibration becomes: 

(8 
y 3 2 x; (A,cos Wt + 8 B, sine t ) (7) 
pes 

In order to calculate a neurmal mode pattern, one cnust resort to 
beam flexure theory. The ample (Zuler) theory of vending may ce assumed 
to hold if the vibration occurs in one of the principal planes of f1¢*ure of 
the car, and if the crossectional di:nensions of the bar are small compared 
to ite length. In the actual case of a ship or submarine the first of these 
restrictions offers no problem, but the second can be unrealistic, particularly 
for a modern submarine with its "cigar"form and essentially circular sections. 


Strietly speaking then, the Luier oending theory is not applicable to a 
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submarine or even a beam whose cross sectional dimensions are 
appreciable. For the preseat, however, we shall use the Euler theory 
as a beginning and modify it subsequently by a correction for shearing 
force and rotatory inertia effects. 

We refer now to Figure A-l, which is a representation of a non- 
uniform beam and an element thereof. The symbols are in accordance 
with the Nomenclature, except that w here is taken as some loading 
in lbsft. With the axis oriented as shown, the well known Euler 


differential equation for the deflection curve is: 


2 

oe ? , 

ety. °° cal - M. (8) 
FIGURE A-1 
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We differentiate (8) once to obtain the shearing force: 


ya 
-——3—— («1-94 ) = - om * -Q, (9) 


An additional differentiation yields the loading, 


eS 2 
3 Oy O@ <« 
~doe— ( 2. 4-4) i w, (10) 
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Equation (10) is the differential equation of a bar subjected to 
a continuous lead of intensity W , Im the case of free lateral vibration, 


ene can appiy D'Alembert's principle and replace W by the ‘inertia load- 





ing 
mm a*y 
Substituting (1!) for w in equatian (16), the general equation for 
the iateral vioration of the var becomes = 
oe ( #1 ay oe dy 
dx” ax yt? 
Finally, since, tor simple harmonic motion,dy_ = - yw 


ot 


we may write equation (12) as follows: 


Oo” (er Q's 
seh Sar) mee 


We now return to the reality that the diinensions of the amas sec- 
tion are not small with respect to the length of the bar, Investigation of 
the effect which the croseé sectional dimensions can have on the frequency 
are incompassed in the classical literature, [6st 9](35]. Of these 
perhaps the most widely accepted is the treatment of the effects of rotatory 
inertia and shearing force by Timeshenko in reference [9 JL 6351. There 
are doubtiess other effects introduced by the cross sectional dimensions 
and configuration but for the purveses of this thesis the Timoshenko 
mechanisw: was adopted because it accounts for tne twe major effects 
mentioned abyuve. 

The effects of beth rotatory inertia and shear deflection become 
of considerable importance in the: higher modes of vioration where the oar 


is subdivided ints shorter segments by an ever increasing number of odes 


[9 jfeo] . 
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a. Rotatory Inertia: 


Refer again to figure A-l. One can easily see how the element 
dx during vibration not only performs a translatory motion, but can rotate 
as well, The variable angie of rotation is equal to the slope of the 
deflection curve, dy/dx, and the corresponding angular velocity and 


angular acceleration will be given by 


_dty_, and _Q3y_ 
Oxdt Oxd.t* 


The moment of the inertia forces of the element about the axis 


through its centroid perpendicular to the x-y plane wiil evidently be: 


dM - 4 (4 ~ a8 dx =\( = m + 2e/4) os =. )( 1, )* (15) 


(The moment is taken positive when it is clockwise), where m/A is the mass 
per unit volume ofthe non-uniform bar, 
We now expand (15), and obtain: 


oy di -_™ * -& @ (m/A) 
(ae) > 9 x 24 2 dx 


£ at . O@M@M/A) 1 \ey (15 a) 


ax (2)* / dx at? 


dM + = = I m/A 


If we neglect second order and higher infinitesimals, 











3 
dM -(-—S =.) dx (15b) 
2 Xd *& 
- 
Remembering . again that for simple harmonic rnotion 2 a = -a)"y, we 
finally obtain: 
mm <--@h 3 Gy) )@ «~ Se 8 (16) 
. A ox A dx 


This moment should be accounted for in considering the variation 
of bending moment along the axis of the bar. Thus, instead of equation (4) 
we will have: 


dM = Q + Ima? dy (17) 
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Sucstituting this in the equation for the deflection curve, 


5S) va 


and using expression (11) and (12), we finally o-taiu: 


{> f. | = ina)" > =F tye Ae oy o*) (19) 


b. Effect af Shearing Force. 


SD still More accurate differential equation ia ottained li not 
eniy rotatory inertia »ut eiso shear deflections is considered, The siope 
of the deflection curve depenaa not only upon the rotation of cress sections 
ef the usar wat on the Gefiection of the sections due to shear as weil, .«t 
“Y Genote the siore of the vac when sheariny foice is neglected and 6 
the angle ef shear at the neutral axis of the same cross section, fvi- 


Gently the total slope must oe the sum of O and, , i.¢@., 


wy 7 7tA [9] 
dx (20) 


vaing the dledeatary vending theory, we Lave for (ending moment 


and shearing force the following expressions, 


a = = peee7 a2 = KACEC = KRfay - “G 
az dx. . 


(2i) & (22) 
K is a nuc erica! factor which is a function of the shape of the 
¢€ross sectional area, usueily defined as fvilaws: 
KR = average spear stress ia section 
manirnuin shear stress in section (23) 
This is the static shear stress distri ution factor used », Timoshenko aud 


discussed in standard texts. 
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The complete mechanism of the bar including both rotatory 
inertia and shear deflection effects may now be formulated. Referring 


to Figure A-l the differential equation of rotation of the element mam n 





Ll 
will be: 
O I asi 
ii 
Ta Res a eee (24) 
Substituting equations (21) and (22) we ootain: 
2 
O KI ge) +e Ge-x) de Im Of = @ 
Ox Ox x A O€? 
(25) 


Neglecting aii but first order infinitesimals, tne differential equation 


for the translatory motion of the same element in the horizontal plane 


will be: 





<z 
— <4 a = oe dx , (26) 


x I< 


ys 
m St - g—bsaa(se a> 0. (27) 


Taken together, equations (25) and (27) compFise the simultaneous 
differential equations which describe the more complete mechanism of 


lateral vibration for non-uniform bars, 


grates: dies polly inal iad/adt in qe elesdpres sae hywes ot 
grivielel Setatacitel ot wna don crete biltedlish wonde bes alieent 
a cit sears oft Neenlienrs Ma oakianpe laitmeyeiibwl b=% rin) Wows 














sed be 
; P tee 7 
(#5) eae is oS . 
va 
— unkesco wm (55) baw (Ih) eemiianpe gelaenireres 
= - , ) cco 
> 
con (ona Ot oe 
ll anal Cr —- | 
“ceutdcupe Jaturesetis 24 eiamentinnink valve tit tot OA gaitnaigalt 
nasig ieoa tod 963 98 sosmmole auae 208 Re mitce retalenane 034 10) 
(eo dae 
a) es ' 
oy eh hn te 
~~» i 


tts) \ 2 _ eet )onal 6. ‘ as nt 


cadassi tenis di diMigma> (Si bus (EX) peeibanpe suilinged wealet 


eer rn 
_ y vad eeebhen cree vu) avibardiy teretal 
— oo w« ail 


—e(= , ,- Es - is Se 


Firnoshenko's aualysis of rotatory inertia and shear deflec- 
tion effects tor a uniform var may be found ia reference [9\ 


He obtains two simulteneous equations similar to (25) and (27). 
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LI ar + &B (2. x) BG - = ot J @s)* 
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He then proceeds to eliminate from equations (28) and (29) v, 


and ovtains thereby the weil known fourth order coupled equation: 
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amit a> , 
KA*G * tr q (30) 


* Timoshenko's notation has been altered to conform with the authors’. 
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Equations (25) and (27) appear not to be susceptizie to such an 
operation, since J, A, m, and X are all functions of x in this case. 
Consideraile effort was spent by the authors in an attempt to eliminate 
Y from equations (25) and (27), without success. If the operation is 
possible, and this the authors doubt, the resulting equation will at best 
be too complex and unwteldly for practical ajplication, 

The difficulties involved in the solution of equation (30) are 
well known. It is quite obvious that the irregular variation of I, m, K, 
and A with X in the case of a ship dictates that equations (25) aad (27) 
tnust be solved in an approximate manner using the graphical, lumped 
parameter, or kayleigh-Rite procedure, The task is so formidacie 
that the authors could find no exact solution in the literature encompes- 
sing both shear and rotatory inertia effects applied to an irregular non- 
uniform) var. 

To recapitulate, the vasic ¢quations used uy the authors are, 
in simplified form: 


f' = - mw’y 


GO = KAGS = BAG ly! - 7) 





m’s Q + mn nT 
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ee7 7 


The Loundary conditions applicable to the quadruy le integration om: 
Q = 9 for = 0, xX & of 
M = 0 fer z= 0, X¥ = 


(34) 


(32) 


(33) 


(34) 


(35) 


(36) 
(37) 
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Equations (31) - (37) descrite an eigenvalue provlem, with Oo being 

the eigenvalue. Solutions wil! exist only for a discrete set of w* . For 
each such value of W* a unique mode pattern will correspond. An excep- 
tion to this is the case when W* = 0, where two possinie modes exist: 
y = 1 (rigid body translation), aad y = 1 - — (rigid Ledy rotation). 
The two modes for this unique case have been orthogonalized for 
convenience . 

Useful ortho,onality relations exist between the normal mode 


patterns. In particuiar, for each 


Lt 
f[ my, de = 9 (38) 
0 3 
ct 
[ =m, ae = 0 (39) 
p 


and approximately, for the lower modes, 


a 


[fm yy@& #0 i; (40) 

v cat 
Murray [71] py application of integration by parts te equations 

(38) and (39) discovered another interesting orthogonality relation. which 


has apparently not been noted previously in the literature, 
=a, x af 
L%:i dx va m (§ - %) d¢ =. 9 (41) 


Equation (41) is used by Murray and the authors to evaluate the constant 
of integration of ft 7‘dx, and this exactly determine the location of the 
true slope axie (This axis has previously peen located arbitrarily at the 


mid-length point of the slope curve). 


* The reader ia referred to Appendix D for the derivation of the exact inter- 
mode orthogonality conditions and for a discussion as to the applicability 


of the approximate condition, Lquation 4°. 
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3 and substitute in equation (41), it is 


may be evaluated as 
© rx 
/ 
“ Wh es. x) dF dx 
3 teas - *®) d§ dx (42) 
0 ~o 
In a like manner, use is made of equation (38) to determine Cy where 
y = + C,. 


It must be emphasized that orthogonality among all the modes must 


/ 
If we take y= N + «@ 


easily seen that C, 


C 


be rigorously snnimecinetl [2] [71 | A starting pattern for the iterative 
cycle is not strictly admissable unless it has been explicitly adjusted to 
conform with the orthogonality conditions described in equations (38), 
(39), and (40). The resultant pattern is not valid until the same adjust- 
ment procedure has been rigorously executed. 

Orthogonality with respect to rigid body translation is insured by 
use of a zero shift of the y-axis: _ * y; + A, Itis clear from 


equation (38) that: Fi L. 


Rf mas s -f ym G&. (43) 
} 0 P 


Orthogonality with respect to rigid body rotation is insured by 
use of a rotation of the y-axis: Y; = i + B (x-x), 


From equations (38) and (39) we may write 
L 


fee Y, j~ - =) a = 0, (44) 


° 
and substitution of the aodove expression in equation (44) yields: 
L 
2 
B je - xk) 2 - |, (x - €) dx. (45) 
f 


Oo 
It is further required that both the starting and resultant deflection 


patterns be orthogonal with respect to each flexural mode prior to the 


one of immediate interest, we formulate the expression 
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¥,2 ¥,+ A + B (x - X) + Cy, + Can teens OEY (46) 


? C, etc are readily determined by substituting %, y, + OM, ) 
etc. into equation (40), For example, if we are interested in adjusting 


Constants C 


the second mode so that it is truly nerrmal to the first mode we set 
ye a, © 5, %, 
whence from (49) 


4 
Cim y?a& = - [» y. y. oe (47) 
0 ‘ ° =! 
In general 
l. L 
, [x ~ ee [e i. de (48) 
i), y, e ; d 


Hot until the overations indicated Ly equation (46) are carried 
out should a cycle be begun nor a resuitant deflection curve se accepted, 
Since the frequeucy depends upon the value of the resultant curve at the 
normalizing point, ¥; *, the frequency showed be calculated only alter 


orthogonalizetion is completed, 
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Part Il 


Detailed Outline of Procedure : 


Thie section supplements the general outline of the procedure 
of this thesis presented in Section I, and sets down the mechanics of 
the operations implementing the relationships derived in the theoretical 
section, Part 1 of this Appendix, 

Seme remarks helpful in interpreting the analytical expressions 
aad the curves of Appendix B are in order, Kornan numerais indicate 


figures in Appendix 4b. 


(i) All curves plotted from data given or calculated are plotted 
at intervals of 0.025 L and faired through ali spots. This is 
considered to give the least number of stations which will 


give a sufficiently accurate determination of true curve shape. 


(2) Since the length coordinate has been ndn-di:mensionalized 

( = x/l), each integration performed with respect to it results 
in a factor of Lin the denominator of the integral expression. 
These factors are rernoved at two points in the cycle: when 
determining curvature, and when normalizing deflection. How- 
ever, their presence within the cycle results in appropriate 
modification of correction ter:ms calculated to add to functions 


in the cycle. 


(3) Since the virtual mass distribution of the starting data was 
im pound units, the m of the theoretical relationships has been 
replaced py mg» the g (acceleration of gravity) has been 
carried throughout as a constant and is removed at the normalis- 
ing stage and in the frequency determination. The g constant 
necessarily apsears also in the various auxiliary functions 


involving m 


= a: 
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(4) The horizontal backstepping of the various integral curves 
in the figures of Appendix B is inherent in the operation of the 
mechanical integraph used. (See Plate 1) For convenience, the 
backstep has been set at integral sumbers of inches, either two, 
three, or four inches, and is readily apparent from the y-axis 
labelled for each integral curve. Coordinates along the horiz- 
ental axis have not been renumbered for each axis shift, >ut 
should be transposed mentally in examining integral curves. 
Vertical shifts of integral curves are made where presentation 


is facilitated thereby. 


(5) Also inherent in the integraph is a scale reduction of the 
integral curve with respect to ite argument, Scale factors used 

in this investigation rauged from thee to 7.5; while not indicated 
directly on the curves, they of course determine the scales marked 
on the integral axes. The continuous maintenance of correct 
scales is vital to accurate determination of frequency from the 


deflection normalizing factor. 


(6) Figures illustrating the method have teen so arranged that 
the four principal integrations comprising a single cycie appear 
in one figure (e.g. Figure VI) while the auxiliary curves per- 
taining to the same cycle appear on one or two separate sheets 
following (e.g. Figure VU). At the end of each mode is a figure 
( XII, XX, XXXI)) illustrating the convergence of the normalized 
deflection patterns from the orthogonalized uniform beam pattern 


ultimately deemed acceptable as the normal mode pattern. 


Sequence of Operations: 


When the tables and srmooth curves (1) of structural information 
have been prepared, it is necessary to obtain the coordinate of the longi- 


tudinal centroid of the appropriate virtuai mass. This may be done by 
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any standard technique, but figure II illustrates a simple and unique 
method. The tangent to the end of the second integral of the maes distri- 
bution is extended to the vase line, and the resulting intersection, rela- 
tive to the origin of the second integral curve, is the abscissa of the 
longitudinal centroid of the mass. The proof of this is given in f7] : 
Sufficient information is now on hand to prepare the tables of 


working functions, including 
mg =% 1/u*a, 10°? x1, KAG/10°, (f- K)/i, 


» 2 = 2 
mg ( & -k)/L. and g fons -R)a g/L 
As they become known, there will be added to the working functions 


table the foliowing tabulated functions: 
¥,™ 4 Yo: Yo 1% , etc, 


This table is prepared for intervals of 0,025 L. 
Certain definite integrals and other constants useful throughout 


the application of the method are also calculated and listed at this point: 


L 
Ayn = [3 dx/ L (49) 
oe 3 
—4e 
N - dx/ lL 50 
. ie x) mg dx/ (50) 
x 
4 
Pp - & ad¢ dx/L 51 
= f[ [lf ® mg ag ax/ (51) 
Subsequently, when data is available, the following ae added: 
-? 
Wis [vag ax/ LL (52) 
“2 
W | Lng . 
5 fy m g dx/L, etc (53) 


It is now possible to commence the process. While any 
approximation to the mode shape can be used ag the starting curve, it 
is found that the most accurate approxi:nation to the final mode shape 


available will lead to accelerated convergence, hence fewer cycles of 
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iteration for the desired accuracy. 

The proposal here is to take a known pattern, as from a similar 
ship or a uniform beam, and force it to compiy with conditions of ortho- 
gonality applicable to the uitimate normal pattern. Thus in the investiga - 
tion conducted for this thesis, the normalized ordinates of the patterns 
found vy word Rayleigh are taken arvitrarily as a starting point. 

The following description traces tne development of a third, or 
four-noded mode, to demonstrate the typical case of a higher mode, 


From the working function tacle and the values of Y_., tabula- 


3 


tions are obtained and plotted of mgt 3° nig (x ~ &) Y 4/ Le, Cagys ¥ 3° and 


mgy,*,- The four resulting curves are integrated (s XII), and the vaiues 
obtained are divided by the operational constauts indicated to ostain 





|e 
A = - Ag Ys toy (54) 
... 
z a 7 L = ~ ~ x) Y {is (55) 
ce 
Cc. # bogs. x i= Me (56) 
- } 
- y ax/f lu 
~. Ama Y, ay (57) 
ii W, 


Now the values of (x - X), J, are multiplied by 5B; the values 


of y, are multiplied by C 


i’ ii’ 


we obtain 


and the values of Y, ere multiplied sy C Then 


ye A+ B + e)/L + Gy + CyYe (58) 


Ali values of fF are divided by y*, its value at ite normalizing point (the 


ecenter-or centermost-after antinode as obtained cy a pilot of that region) 
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to obtain a tabulation of y , the orthogonalized, normalized starting 
curve for the first eycle of iteration. (But see discussion of orthogo- 
nality in Appendix D). 

The cycle itself (KX HI) commences with the curve of inertia 
loading = - ymoQ » which is calculated, plotted, and integrated, to 


give a function of shear, 
gab ~y¥mg (59) 


where the constant of integration is zero, since the integral curve 
starts at the origin, and the boundary condition of zero end shear dictates 
that the true shear curve should also paas through the origin. However, 
the right end of the integral curve may differ from sero by a small 
amount (usually iess than 5% of the maximum ordinate of the curve). 
When this occurs, a linear correction is applied to the curve to again 
conferm to the above boundary condition. This is done by drawing a new 
axis joining the ends of the shear curve, measuring its ordinates from 
this curve but norma) to the old axis, then replotting on the old axis. 
This correction is illustrated in Figure Vi, where the reader should 
compare the curves marked gQ/\ and g0/L (corrected), but has been 
omided from the presentation elsewhere to avoid congesting the plots un- 
necessarily. In the first cycle, the corrected curve gives gM'/i =9Q/ lL. 
(At this point in all but the first cycle, each value of shear receives 
an additive correction to account for rotatory inertia, the correction 
being a function of the bending slege 7 obtained from the previous cycle 
in a manner to be described, This slope, in the form L?Y , is multiplied 
by the working function mg 1, A and the product added to the shear 


term to give a curve of 


gM /L = gQ/i + Ye m1, LA (60) 
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This is the only procedurai respect in which subsequent cycies differ 
frorn the first). 

The curve of 9 M'/\L is integrated, and the boundary condition 
ef sero end moment is applied im exactly the manaer described avove 
for the shear correction. In both cases the correction can be considered 
to filter out accumulated errors resulting from smail inaccuracies of 
integraph operation and calculation. The values of & i. ad a direct 
function of moment in the hull girder, are tabulated from the corrected 
curve, and each is multiplied by the corresponding value of 10°, kl, 
to give a function of curvature, g7nlo’, which is plotted. 

Curvature is integrated to give bending slope to within an 
additive constant. To this component is added the shear component of 


siope, 
gf x 10° L = 3 Q/ i) + (%AG/10°) (61) 


based on the shear vaiues ootained after correction of the first integral 
in the cycle. 
The addition gives a function of total slope, still without the 


constant of integration: 
’ 6 —~ 6 F 2 & 
cy x10 /L ® Ze x 10 /L + gr i@ si (62) 


To determine the constant, or axis shift, use is made of Murray's 
formulation (see Part I of this appendix) on an auxiliary plot (XXIV). 
The values of 9 ’ 10°, lL, from the above addition are multiplied by the 
working function 


. . 
[8 ee > BMce 7h 


and the product plotted and integrated on the auxiliary plot. The value 
of the definite integral obtained, ) 


6 fonino® feng (+ xd dx/L 
7 3 5,/> = <« : ee - £ ous (63) 
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is added to the uncorrected total slope to give the function of true slope: 
F 6 f 6 46 
By xl0s;L =# Q7 ~lO/lL + Wg C4/b (64) 


This is plotted on the cyvle plot and integrated, yielding, to 
within an additive constant of integration, a function of deflection, 
SJ yx 10°, nu To obtain the constant, this is multiplied by ™g 
and the product plotted on the auxiliary plot (XXIV) and integrated. 
The value of the definite integral is divided ty a function of displace- 


ment to give the constant 


10°g Cyn es ao t gers ie) me ax/L (65) 

This is added to the values of gyn 10°, ‘, to give the 
deflection pattern 

Qy * 10°, x? = 9% 10°, hy + 9 Cy x 10°" (66) 

This pattera is orthogonalized in the sarme manner as was des- 
¢ribed for the Kayleigh pattern, to give gy* 10°, hs The value at the 
normalizing point, ay" ie / ”. is recorded, and divided into all 
ordinate values, to yield the normalized deflection pattern, Y 3-4 ; 
of the first cycle. it is this pattern which is used to compute the inertia 
loading for the next cycle of the iteration. 

The value of BY - « 10°, 4 is multiplied by i”, and the 
product resulting ia dividedinto the values of QO 'x 10° from the last 
cycie of the iteration for each snode to yield the normalized curvature 
for the mode, which now wiil have the correct relationship to the 
normalized deflection pattern. 

The value of ¥ corresponding to the cycle just completed is 
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The value of W ° cornputed from the last cycie performed is accented 

as the natural frequency corresponding to the mode being considered, 
In preparation for the rotatory inertia correction, the shear 

curvature is suotracted from the total curvature to give the corrected 


bending curvature 


gx x10") = 97’ x 10° / 4. - 98 x wesL (68) 


which is aormalized by dividing by ay > = 10°, »* to yield the 


required vaiue of J L.. 
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APPENDIX EF - SUMMARY OF DATA AND 
CALCULATIONS 


Figures I through XXXIV in this section depict 
the application of the method developed in this thesis 
to the model for which characteristics are tabulated 
in Appendix k.. The calculation of argument curves 
plotted fron: data is iliustrated in the sample calcula- 
tions of Appendix C. 

Table Il summarizes the effect of each iterative 
cycle, and its associated orthogonalizatios check, oa 
the natural frequency of the mode considered. The 
corresponding effect on deflection patterns is presented 
in Figures XU, XX, and XXXII, for the first, second, 


and third modes respectively, 
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APPENDIX ¢ 


Sarmple Calculation 

This appendix contains a tabulation of the working ‘unctions for 
iteration, a table of constants used in the iteration, and a sample 
calculation in tavular form, AIL symbole are in accordance with those 
appearing in the Nomenclature, 

The working functions (Table IV) pertain to all of the modal calcu- 
lations, These functions were calculated asa matter of convenience as 
wellas necessity. their form is governed both by theory and wy the 
conventions adopted and explained in ..ppendix A. 

Table V contains a list of constants which, in general, apply to 
all three medai calculations. Orthogonality adjustment constants W, and 
W, apply to all meda! calculations beyond the first mode and the second 
mode respectively. ‘Ihe derivation and explanation of the constants arising 
from integration may be found in Par! 1 of Appendix 1, The usual value 
of Young's Modulus was adopted, and the shear modulus was then computed 
from the well known relation: 

. ee cr ae 
where wp is Poisson's hatio, taken equal to 0.3. 

The sample calculation embraces the third cycle of the third mode, 
This mode was chosen as cteing the most generz! of the three with respect 
to the orthogonality conditions, Ihe third cycle of this mode was chosen 
for the calculations because it includes the rotatory inertia correction, 
(The second cycle does also, but the results thereof! are less conclusive 


than those of the third cycle). The calculation is presented in tadluar form 


in Table VI. The source of each function appears beneath its symbol at 
the head of the column, Amplifying remarks pertinent to the sample 


calculation appear at the end of Tabie VI. 
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CONSTANTS KEQUIKRED FOR ITERATION 


CONST ANT SYMBOL VALUE 
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i Gx 
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[ [Ez 6 ate 
Oov6 Le hes Le 
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[oz mg dx N 399 1b/ft 
t) L* L, 
"2 
[3 mg dx Wy $062.5 ft-lb 
. is 
‘. : 
fo mg dx W 5040.9 ft-lb 
a 4 a ‘ 
Length of madel L 160.5 ft 
Longitudinal centroid ii (0.472) L ft. 


(from bow) 
Shear Modulus G 1.5 x 10° te /in. * 


Young's Modulus kL 30.0 x 10° tejte.* 
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TABLE VI (Continued) 


EMARKS RELATIVES 19 SAMPLE CALCULATION 





COLUMN REMARKS 

3 <--- e-ereee-- Inertia loading with cd* taken equal to 1, 

ene wn ewcewenne %.248 the normalized component of slope due to 
bending, taken from the previous cycle. 

ees 3-2 3*ML is the derivative with respect to x of 

LA 
the moment of the inertia forces, st any section, 
about the vertical ceatroidal axis. When G2 3) ms 
Lon 

is added to the shear force, one obtains dM/dx, 
which re‘iects the effect of lateral inertia forces 
plus shear due to rotatory inertia within the 
approximation imposed by using YW from the previous 
cycle, 

Q anew cee n anne -~-+ Component of slope due to bending moment deiore 
determination of its true axis. 

IU -------------- Component of slope due to shear flexibility. 

Ih ---+------- ---- Column 9 plus column iG. Slope before true axis 
determination. 

6/*% ” 
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TABLE VI Continued 


COLUMN REMARKS 


“6 
15 10 4 gom dx = ~- S175 
f a 7 





t 1° m 
¢ (ase) ° aE = 
+ S- a 
-5175 
az -« “$45).~—COS 4 ¢ 0. 6023 
17, 18, 19, 26 Calculations necessary for determination of 


adjustinent constants to be applied to } °3,Ya3 
in order that the resultant deflection patiern is 
orthogeoqal with respect to: 
1, Rigid bedy translation 
2. Rigid body rotation 
3. First mode deflection pattern 
4, Secoad mode deflection 
pattern 


17 i. Urthogonality with respect to rigid vody translation: 
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Condition: [yume = @ 
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JABLE VI Coatinued 


COLUMN ZMAKKS 


18 2. Orthogouality with respect to rigid body rotation 


& , 
Condition: [ume (x -x) dx = 0 
d 


‘Ne@eault: 


a 
6,- ~ 
10. 2 3-3. x - X) dx = -44.8 
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L - . 
10°82 Ya. m (X-*) dy 
i? 


“Biu) = - £ 


N 


-44,8 = + 0.1493 
300 


19 3. Orthogonality with respect to lat mede: 
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Condition: ly, y,™m ¢x = 0 
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29 4, Orthogonality with respect to 2nd mode 


L 
Coadition: Pot y, m dx «= 0 
6 


Sot 


bawelenss £Y Rumer 


amaaaa 


4 ae-se ‘ge b2ilban 
eet 


wats tbe Eo a 


2 Gs) ok yor » om 


oe ae 


— uieot Tel ef foaqace Gilw 
bert yilawegeaert 4 


~——— 


— a) quam (ms aner, ees 











403 


TABLE Vi Continued 





COLUMN REMARKS 
“62 . 
29 Result: Og y ym @ = -1l2 
ee == L 
L 
Jog? F m dx 
LEI 
. a W 
"2 
-112 = + 0.0222 
5940 
. on | Bae 5 
2\ Adjustment to ve applied to 10 g ¥;.; ior 
og 
orthogonality with first mode deflection pattern, 
6 
22 Adjustment to be applied to 13 y ier 
3 PP ——3 ae 


orthogonality with second :nede deflection pattern 


6 
23 Adjustment to be applied to ss i Bs for 


orthogonality with rigid sody rotation. 


6 . on .. 
Q , . ae 
24 8," A etyy! A + B(x-x) + CY; + C., Ys 


The value of 10°, % at the positive anti-node 
— © 


ASA (near X/1. * 0.5) 


ia the normalizing value aad also a function of the 


eigenvalue, Its value is 3.1503 for this cycle, i.e., 


6 
log y* = 6.1503 
we oe 
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TABLE VI Centinued 
REMARKS 


Norrnalized resultant deflection pattern, Y5 . 
, “7” 
1G . 


330 i (0.1503) 


Normalization of Lo for use in next cycle, 


CALCULATION OF FKEQUENCY 


6 2 6 


l . 20 g/i _ ie” 50.8 
RL. aL NS eg ote: 7 2.14 a. 68 sk 
“ iG” gy* 3.1503)(100. 5 
a We = ( y ) 
3168, 025 
"1503 : oe 


145.64 rad sec 


2™ . 6, 28318 x 1000 = 43.14 Milli-seconds 


35.5 145.04 


60 = 145,64°69 = 1391 c.p.m, 
"RR &, 2838 
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APPENDIA OD - SUPPLEMENTARY DISCUSSION 


l, Secondary Lifects Not Treated: 


Shear Lag: The classical beam vibrations as applied to ships 
do uot take special notice of the fact that resistance to deflection is 
concentrated generally at the extremes of the cross section and that, 
in horizontal bending, stress transmitted from the deck to the side will 
not maintain a constant leve! down the side as indicated by the constant 
distance from the neutral axis. This is because of the effect known as 
shear lag. While important in lightly framed ships where decks and 
longitudinal culkheads play little part structurally, it was not thought 
to be significant in the ring-frarmed model of this investigation. Shear 
lag is treated in detail by K. A. Anderson and J. C. Houvolt ("The 
Effect of Shear Lag in the bending Vibrations of fox Beams", NACA 
TIN No. 1583, bay 1948) and others. 


Second Frequency Spectrum: The simple Harmonic assumption 
applied to Timoshenko's complete mechanis.n leads to an equation which 
ie fourth order in circular fequency , indicating the possibility of a 
second reali value of frequency for a given mode, Traili-Nash and 
Collar [63 | were able to determine this second spectrum of frequencies 
for certain special cases of uniform beams. Its physical importance was 
not determined, but it; was thought to arise from a resonant interaction 
between rotatory inertia forces and shear stiffness. 

The iterative method of this paper, while considering both 
effects, and permitting them to interact in the sequence, is based on a 
singie assured frequency, and does not facilitate simuitanecous inter- 


action. Hence the second frequency spectrum is not explicitly noted. 


Torsional Coupling: It should be noted that coupling can occur 
between torsional aad horizontal flexural vibration, with resultant effects 
on the patterns and frequencies of each, Such coupling has aot been treated 


in this investigation. 
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2. More Exact Orthogonality Formulation: 


A review of the theoretical development of this method revealed 
that the formulation of the orthogonality relationships used were inexact 
when applied to orthogonalizing with non-trivial modes when rotatory 
inertia is considered. The effect of the discrepancy is discussed in the 
next sucsection. The development of the more exact formulation® 
will now be given. 

The normalized deflection and bending slope, and the circular 


frequency of each of two different modes are given vy 


¥y Bs G) y; 1% 3 GQ); - 
Multiplying equations (25) and (27) of Aspendix A-1 vy uF and y; 


respectively, and integrating, we have ne pee t : 


0 


L 
t Im <2 ; / 
{ (21 ¥; idx + (4 YIidx + [xa Y; iF & J. j )jdx = 0 (69a) 
= o 


re 
L , P 
fos Y;¥j,Pi ae + [waG (y, -%;)] ¥, @& = 0 (69%) 
0 o 
and using the multipliers from mode i, we have for mode j , 
tu. u E 
2 





/ 
LY. im W. ¥.8; 'y, | 
face )dx + 5 3 ee + fanary Yi - 3; Be Yes «= 9 (70a) 


0 o 

L 7 : * 

m Y, ¥; ff dx + [RAG y; - yy )| Y. dx = 9g (70d) 
(4) 0 


Integrating the first term of (69a) sy parts, we have 
Lc 


LL ls 

. (ee im ,< oy ! 

el Y, 5; | . ic ¥- uo dx + [eta + KAG(y/y -¥.¥-) 
; i aor 0 (71) 


*DOue to Prof. 5. H. Crandall, Massachusetts Institute of Technology 
andto ©. u, Dolph, Willow Num Aesearch Center, University of Michigan 
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where the first term vanishes since curvature in a free-free beam 
vanishes at both ends. The corresponding operation ie carried out on 


(70a), and the result sudtracted from (71), giving 
L ib 
2 3 Im ' ' 
(wh ob) fa 3 dx + [ wae ly: v - ¥; % | dx = 9 (72) 
Now (69) is rewritten, integrating the second term b, parts, giving 


fo y. iF a? dx + KAG Cy: - v; yi} - [easty'y -% yj; eda = @ (73) 
where the second term vanishes, since KAG (y. - 7. ) is identical to 
Shear, which must vanish at both ends of a free-free beam, The corres- 
ponding operation is carried out for (705), and the resulting equation is 
subtracted from {73), giving 
( Ww." - ws") f my: ¥, dx + [se | 7% ye +e ¥: | dx = 9 (74) 
Adding (72) and (74) gives the orthogonality relationship 
© Im 
(ar 2) [ (4% + my, x, ) ax = 0 (75) 


which, for ww, # ud: » leaves the useful relationship 


, 


[x7 + my, St )dux = 9 (76) 


This should be compared with the relationship used in the investigation 


[vey ¥, ax “— (77) 
3 Z 


which is seen to be in error by the integral of the first term in (76), 

The importance of this term varies fram mode to mode, and 
must be evaluated for higher modes to determine the necessity of 
applying the exact condition, It does not lend itself readily to application 


in the iterative procedure of this thesis, since it provides corrections 
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to bending slope and deflection, quantities which vecorme available at 
different stages of the process, and which cannot Le corrected sesara- 
tely. Aa acceptable, though somewhat unwieldy apsroach would ce to 
apply the complete check at the deflection level, examining the slope 
correction te see if its magnitude warrants its use. If it is significant, 
the Lending slope is corrected, and the new total slope integrated. The 
new deflection curve is given the exact orthogonality check, this time 
assuming that the sle,e needs no correction. 


In the case of this thesis, an auxiliary investigation showed 


that the complete check changed the natural frequency for the third modes 


(previously corrected b, the approximate relationship) by 0. 36%, and 
that the effect was even smaller on the second mode frequency. It 
appears that for the modes ofinterest, the added inconvenience of the 
full correction may nave but little justification. Where applied, the 


correction is stated vy the following equations: 
(Y5% r yj _.- =e (78) 


(F; ), = iF + CY, + Coe “ (79) 


where the constant for each prior mode is defined by: 
L ” L 
{= 0 Odx + om 5 dx 


. - mL ox sae y , dx (@0} 


é 


3. Third Mode Frequency: 
As noted in the Discussion of Results, the natural frequency 
obtained for the third mode (1391 cpm) oy the iterative process and 


orthogoaality flliering is considered to be slightly high. This section 
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presents estimates of the correct value ontained from several independent 
approaches, 

a. The cycle-ty-cycle behavior of frequencies in the third made 
caiculation manifestsalternation casting some douct on the value from the 
third cycie. Consideration of this behavior seems to indicate ultimate 
convergence to a value somewhere between the second and third cycle 
results, whereas carlier modes showed continually falling trends. The 
change in sehavior in the third mode appears to stem from a starting asa- 
sumption which was in error, even after orthogenalization, to a greater 
degree than were those of the other modes. This may be due in vart 
to the use of the inexact orthogonality check, as well as to any of the sources 
of error listed in the Discussion of Results. We conclude from this approach 
that the third mode natural frequency should be in the neightorhood of 
1385 to 1396 cps. 

b, The nore exact orthogonality check formulated in equation 
(80) was applied to the result of the last iwration of the third mode, dis- 


regarding the effect on slope. Tne constants Secame 


-. 605 - 27 | 
© * - Seep sas7.8 * * %-eb4es 
instead of 
- 595 
"es rrr =z + J, 0626 


and 


-1iZ + i 2 oy 
e : : + 3.92166 
Ci ’ 540 + 76.3 . 


instead oi 
-112 


i = ~ “Shan + 0.3222 


oC 


The net effect was to lower the natural frequency from 1391 to 1366 cpm. 


The change in mode shape wae insignificant. 
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ec, The empirical ratio of natural frequeucies observed for a 
numberd@ ships, 1: 2: 3, although approximate, and not verified for 
submarines, is applied for interest, It results in no significant change 
in the second mode, but depresses the third mode frequency to 125] 
cpm. This is clearly much lower thaa there is reason to expect, cut 
tends to verify the suspicion that 1391 is perheos high. 

ad. Table VII lists frequencies for the model under considera - 
tion, obtained theoretically. Tabulated with the values ostained formaliy 
in this investigation (MIT) are those ovtained by David Tayler Mode} 
Basin (DTM) by analog computer and those obtained by the authors at 
Underwater Explosion Kesearch Division (UL KD), both sets including 
shear flexibility but not rotatory inertia. 


TADLE VII 
CALCULATED MODEL FREQUENCIES (cpm) 
ode MIT DTMS VEKD 
I 417 444 42% 
é &38 888 &57 
3 1391 1392 - 


Since the added refinement of rotatory inertia should tend to lower the 
frequencies of this thesis as compared te the other sets, it again appears 
that 1391 cpm is too high for the third mode, and that perhaps a reasonable 
value to expect would be 132) cpm. This is based on both the ratio series 
for the DT Mii data and the suspected trend of equal or greater reductioa 


with mode. aumber due to rotatory inertia. 


I¢ shoud be pointed out that neither set of comparison frequencies 
is corrotorated, and experimental frequencies have not yet been ovtained 
for the model under consideration. 

e. Considering all of the above factors, it cau best be concluded 


that the third mode frequency is 1355 + 2.6%. 
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APPENDIX E ~ STRUCTURAL DATA 


Table VIII contains a tabulation of structural data necessary for 
the modal calculations. The source of these data is the Underwater 
Explosions Research Divison of the Norfolk Naval Shipyard, under whose 
auspices the model submarine referred to in this work was constructed. 

The mass appearing in Table VIII] includes the structural mass 
of the submarine, the mass of her ballast water, the mass of her solid 
contents, and the virtual water mass. 

The virtual water mass was calculated from the following 
formula: 


my, bias * (81) 


where b is half the depth of the submarine structure at any section, ' 
including the protuberance of the rudder and conning tower where appli- 
cable, and e = density of water in lb/ft ¥ 

Equation (81) is based on the assumption that the cross section 
is elliptical with the major axis congruent with the vertical center line 
of the submarine, The wash of the water through the superstructure is 
neglected, The resulting curve of virtual water mass was further modi- 
fied by an empirical correction explained in David Taylor Model Basin 
Report # 632. This correction consisted of multiplying the m, curve 
by 0.83 over most of its length, with a more severe correction at the 
ends. Only one-third of the additional virtual water mass contributed 
by the rudder and the conning tower was used. 

These corrections, a Somewhat arbitrary, as all estimates of 
virtual water mass tend to be. The total empirical correction causes 
only a very slight change in the ship vibration behavior. No experimen- 
tal data that verifies the virtual water mass distribution is available 
for this submarine. 


The ends of the total mass curve used by the authors were altered 
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slightly so that the ordinates did not approach cero rapidly, out remained 
at an arbitrary small amount, dropping abruptly te zero at X= 0 and 

X= L. This was done to conform with a similar adjustment to the rnom- 
ent of inertia curve. The need for the corrections arose when, in their 
absence, spurious values of curvature were obtained at the ends Lecause, 
in general, the curve of inertia tended to converge on zero more rapidly 
than the curve of bending mement. 


The K factor used to determine 6 from Q is defined as 


x « average shearing stress in cross section 


maximum shearing stress in cross section 


al 1% 


** _—- (32) 

where I] and A have already been defined, and 5B is the total thickness 
of steel through which the vertical center line passes. (The center line 
vertical and flat keels were exciuded frorn the calculation of K). 

The integral of equation (42) was caiculated separately for each 
component of the structure (i.e., the pressure hull, outer hull, deck, 
etc. ), and all the integrals sumrned to give a value applicable to the 
entire section. Only members which inake major contributions to the 
submerine's longitudinal strength were considered in the adove calcu- 
lation aad in computing sectional moments of inertia. A discussion as 


to the validity of this practice appears in Appendix D, 
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